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We study two double dot systems, one with dots in parallel and one with dots in series, and 
argue they admit an exact solution via the Bethe ansatz. In the case of parallel dots we exploit 
the exact solution to extract the behavior of the linear response conductance. The linear response 
conductance of the parallel dot system possesses multiple Kondo effects, including a Kondo effect 
enhanced by a nonpertubative antiferromagnetic RKKY interaction, has conductance zeros in the 
mixed valence regime, and obeys a non-trivial form of the Friedel sum rule. 
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The continual advance in the ability to engineer de- 
vices on the nanoscale level has led to the recent fabri- 
cation of double quantum dot (DQD) devices [l|, [3, H, 
l1^uIiJ3L both from semiconducting heterostructures 
JJ, |J,HnJHj @| and from carbon nanotubes 0, Si ■ Like 
their single dot cousins [^, , these devices are highly 
tunable. This tunability makes DQDs both a leading 
candidate for a solid state realization of a quantum qubit 
and an ideal laboratory to observe strongly correlated 
Kondo phenomena. 

Single dot devices have long been known to exhibit 
Kondo physics 0, By the use of a gate voltage, 

the number of electrons on the dot can be adjusted so 
that it is odd. By virtue of Kramers degeneracy, the dot, 
through hybridization with electrons in connecting leads, 
becomes a realization of a single impurity Kondo system. 
Transport measurements are the most striking signatures 
of Kondo physics in such devices, where, for example, the 
observed finite temperature linear response conductance 
matches that predicted theoretically 
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With different possible dot geometries, Kondo physics 
in DQD devices is necessarily more rich. Dots arranged 
in series have been shown to exhibit Kondo physics in 
competition with an effective Rudcrman-Kittel-Kasiiya- 
Yosida (RKKY) coupling between dots 0, i, [H, Q [ll] . 
And in parallel dot systems, Kondo physics has been ob- 
served 0,01 where it is expected to appear in conjunction 
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with interference effects 

In this letter, we present a powerful theoretical ap- 
proach to DQD systems. We study two models of such 
systems, one with dots in parallel, one with dots in sc- 
ries. We argue both models admit exact solutions via the 
Bcthc ansatz. The model of dots in series possesses an 
SU(4) symmetry [l3| and is relevant to the study of the 
Kondo effect in carbon nanotubes [l^. For the model 
of dots in parallel, we explicitly demonstrate that trans- 
port properties may be extracted from the exact solu- 
tion. Transport in the system of parallel dots exhibits 
two interesting features: i) at the particle-hole symmet- 
ric point, we find a cooperative combination of Kondo 
and RKKY physics where a non-perturbative antiferro- 



magnetic RKKY effect serves to mediate the formation 
of a Kondo like singlet; and ii) in the dots' mixed va- 
lence regime, we discover a form of the Friedel sum rule 
where contributions to the scattering phase come both 
from electrons in the dots and the leads. 
Models Defined: We examine two generalized Ander- 
son models coupling two leads {I = 1,2) to two DQDs 
(a = 1, 2) in two distinct configurations. In the first 
configuration (denoted by PD), the dots are arranged in 
parallel. The corresponding Hamiltonian is taken to be 



Til = Tin + 



Hunt with Ho 
and 
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Hunt = ^ Vlaiclidcra + h.C.) -|- ^ Uaa'n^ania- (1) 
l(Ta aa' 

Here the cia/daa specify electrons living in the leads and 
the dots. Via measures the tunneling strength between 
the dot a and lead I. Uaa' characterizes the Coulom- 
bic repulsion between electrons of opposite spin living on 
dots a and a' . 
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FIG. 1: A schematic of two DQD systems. 

In the second model (denoted SD) the dots are ar- 
ranged in series. The corresponding Hamiltonian has a 
similar form to the first H2 = P{Ha + H2int)P- Here 
H2int is of the form 

H2int = ^ Vi(c3'^(icr( + h.C.) -I- ^ riaana'a' 
(7,1 [aa)^{a' 17') 
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P is a projection operator which forbids total occupancy 
of both dots from exceeding two electrons. In a regime 
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where fewer than two electrons sit on the dots, the role 
played by both the projection operator and the corre- 
lated hopping terms will be minimal. We note that 
the addition of an interdot hopping term, i.e. (57^2 = 
t±{d\d2 + dldi), does not spoil integrability 21 1. 
Integrability of PD: To analyze PD (two dots in par- 
allel) we first map the problem to an Anderson model 
involving a single effective lead. To do so we need to 
assume the ratio of left /right lead couplings are equal, 

i.e. VLa/VRa = VLa'/VRa'- Writing Ce/o = {VL/RaCL ± 

Vr,lc,cr)/V^, with r„ = {Vl^ + the Hamil- 

tonian factorizes into an even and an odd sector, Tii = 
Ti-e + Ti-o- Only He couples to the dot and is given by 

T-ie = HeO + T-ieint with HeO = ~i I^oo dxcl^d^Cea + 

J2aa ^da-riaa and 
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In contrast, Tio is trivial: Tio — Heoi^ ^ o,tda — 0). 
The integrability of PD is then equivalent to the exact 
solvability of Tie ■ 

To determine under what conditions Tie admits eigen- 
functions of the Bethe ansatz form, we begin by comput- 
ing both the one and two electron eigenstates. The one 
particle wave function appears as 
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Solving the Schrodinger equation, Tie IV') = ^IV")? 
wc find ga{x) ^ 6'(x)e*«^+*^/2 + 0{-x)e'''='-'^/'^ with 
6{q), the impurity scattering phase, to be d{q) ~ 
— 2 tan^^(^^(rQ,/(g — e^a)). To compute the effective 
scattering between electrons, wc study the two particle 
eigenfunction with spin projection, Sz ~ 0: 
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Again solving the Schrodinger equation Tie IV') = 
(g -I- p)|V') gives gixi,x2) to be g{xi,X2) = 
gqixi)gp{x2)4>{xi2) + ixi ^ X2), witha;i2 = xi-X2. Here 
gq/p(x) are one particle wavefunctions with energies q/p. 
4>{x) governs the scattering when two electrons are in- 
terchanged. It takes the form 4){x) = 1 -I- i7(g,p)sign(a;). 
We find that ^{q^p) is consistently determined to be 



7(9, P) = 
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where c and c' are a-independent constants or 



Tn 



(8) 



The first set of conditions describes single level dots ab- 
sent an interdot coupling while the second dots with de- 
generate levels but with a highly finely tuned interdot 
interaction. We will thus focus on the first. 

Exact solvability is predicated on how 7(9, p) deter- 
mines the scattering matrix of the two electrons. The 
scattering matrix has the general spin (SU(2)) invariant 



form, = bip,q)I^l'>' +cip,q)P-,'>\ where a, 6 =T,i 



and I/P are the identity /permutation matrices. The co- 
efficients, b{p,q) and c(j),q), are determined by j{p,q) 
from the relation, q) ~ c{p, q) = (j){x > 0)/ V)(x < 0) = 
i7(p, g))/(l-i7(p, q)), together with q)+c{p, q) = 
1 which arises from considering the eigenfunction in 
5z = ±1 sector where interactions are absent. 

In order for the Hamiltonian, Tie, to be integrable, a 
minimal condition is that the above S-matrix satisfy the 
Yang-Baxter relation. The Yang-Baxter relation governs 
the scattering of three electrons and it enforces the equiv- 
alency of different scattering orders. It is well known 
in the case of an SU{2) symmetry that the validity of 
the Yang-Baxter relations is equivalent to the condition, 
b{P:l)/c{p,q) - 9(9)), where g{p) is an arbi- 

trary function |22l |. In the case at hand, g{p) is given by 
g{p) = {p- - [/„/2)V(2r„C/„). 

Having determined under what conditions the dot-lead 
Hamiltonian is exactly solvable, we are now in a position 
to construct A^-particle eigenstates in a controlled fash- 
ion. An eigenfunction with spin, Sz ~ N — 2M, is char- 
acterized by a sea of N electrons each carrying momenta 
{Qi}iLi and so total energy E = J^i Qi- Iii a periodic sys- 
tem of length L, integrability allows us to write down in a 
compact form the (/i-quantization conditions (the Bethe 
ansatz equations): 



JqjL+iS{qj) ^ 
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These equations are identical to those for the ordinary 
Anderson model but for the form of d{q). The M 
Aq's appearing in the above equations are indicative of 
the spin degrees of freedom. 

There are two integrable generalizations of the PD 
model: i) N-dots in parallel and ii) a pair of dots in a 
T-junction. These systems, to be integrable, must sat- 
isfy a set of constraints similar to Eqns. 7 and 8. 
Integrability of SD: We demonstrate the generic inte- 
grability of SD along similar lines - it was already known 
to be exactly solvable for infinite U [i^. Here, however, 
wc do not transform to an even/odd sector. Wc thus have 
four different types of fcrmions/dot degrees of freedom. 
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This will ultimately lead to the model having an exact 
SU{i) symmetry at all energy scales. 

Constructing the two particle eigenfunctions as before 
leads to an 5f7(4) S-matrix provided Vi = V for all I 
and eda = for all a. (This latter constraint can be 
relaxed enabling one to study non-degenerate dots where 
the SU{4:) symmetry is broken, say by a magnetic field 
or a gate voltage.) The S-matrix takes the same form as 
Eq. [7] but with a, 6 one of four values, (1, t), (1, j), (2, t 
), (2,i). Again we have b{p,q)/c{p,q) = i{g{p) - g{q)), 
sufficient for the Yang-Baxter relation to be satisfied, but 
with g{p) = {p-ea- U/2f/{V^U). 

To construct the N-particle eigenfunctions, we again 
employ a nested Bethe ansatz jl^] . Crucial to these eigen- 
functions being of the Bethe form are both the projec- 
tors, P, and the correlated hopping term (Eq. [J) of 7^2 ■ 
This demonstrates that in the case of dot systems with 
orbital degeneracies and a finite U Coulomb repulsion, 
and unlike Hubbard models with orbital degeneracies, it 
is possible to find a simple Hamiltonian which is exactly 
solvable. As was demonstrated in Ref. ji^l, the {N > 2)- 
particle wavefunctions of any simple finite U orbitally 
degenerate Hubbard model are not of the Bethe type. 

The quantization conditions of an N-clcctron state car- 
rying momenta {qi}fLi arc of the form 



fc-i \ fe \ fe — 1 I i 
Aq, — +2 



n 
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Here ^2(9) = -2tan"i(l^V2(g - fd))- The quantum 
numbers, {A^}, fc = 0, 1, 2, 3 (with A" = gil)) correspond 
to both spin and orbital deg rees of freedom. In a prelim- 
inary analysis of Eqns. [2l|, we have verified the SU(4) 
Kondo physics expected in a regime where one electron 
sits on the dots 18|. Thus the marginal correlated hop- 
ping term in Eqn. Q does not influence the universality 
class into which the physics falls. We also note that the 
integrability of SD can be generalized to N-dots arranged 
in triangles (N=3), squares (N=4), etc. 
T=0 Conductance of PD: In the remainder of the 
paper, we focus on extracting the features of the T ~ 
linear response conductance of the parallel dots (PD). 
The structure of the linear response conductance, G, for 
the parallel dots is much richer than that of a single dot 
containing distinct Kondo effects, novel applications of 
the Friedel sum rule, and quantum critical behavior. 

To compute the linear response conductance of the 
dots we closely follow Ref. . The approach is based on 
the observation [20| that the impurity scattering phase, 
<5imp, of an electron is determined by the shift in the elec- 
tron's momentum due to the presence of the impurity, i.e 
p p + Simp/ L. From the Bethe ansatz the full momen- 
tum of the excitations are readily extracted. Then isolat- 



ing the term in the momentum scaling as the inverse sys- 
tem size, L~^, allows the Jimp to be computed. In the ge- 
ometry we have chosen, G is given by 2e^//isin^((5ii„p/2). 
For a detailed description in the context of PDs see Ref. 



2 

1.75 
1.5 

1.25 

0.75 
0.5 
0.25 


-0.25. 



1 1 1 1 1 1 1 

■ U| = l.U, = 0.5 


t, 

'\ Standard 
1 1 Kondo 


. r|=o.o5Uj 


1 1 displaced electrons 


" r,=o.iu, 


1 I - - G(e'/h) " 


- Ae = 0.25U, 

■ 




■ ^ 






\ > ^ : 


■ ; 

■ 1 RKKY ' 
1 Kondo 


V ^ ^' \ 


non-trivial "^^^ J 
Friedel sum rul^^ 
... 1 ... 1 ... 1 ... 1 ... 1 .. . 



RKKY Kondo 

lix 



Standard Kondo 



FIG. 2: The total conductance and number of displaced elec- 
trons per spin species vs. e^i , Ae fixed, for two dots in parallel. 

We plot an example of the linear response conductance 
for two dots in parallel in Fig. 2. We consider the asym- 
metric case where Ae = e^i — 6^2 ^ ri,r2. The con- 
ductance is plotted as a function of e^i keeping Ae fixed 
from the particle-hole symmetric (p.h.s.) point of the 
system (i.e. C/q, = —2eda) where two electrons sit on the 
dots to a point where both dot levels are well above the 
Fermi level and the dots are nearly empty. In this plot, 
the structure of the linear response conductance, in com- 
parison with that for a single level dot, is more complex. 
This reflects both the presence of interference as well 
as distinct types of Kondo physics. 

At the p.h.s point, (i.e. edi/Fi = —10 in Fig. 2) 
the two electrons residing on the dots together with elec- 
trons in the leads form a singlet. With no bare di- 
rect exchange, singlet formation is mediated solely by 
virtual hopping processes which here promote antiferro- 
magnetic correlations. We thus term this formation the 
RKKY-Kondo effect to mark the role of electron itiner- 
ancy. This phenomena is distinct from the ferromagnetic 
RKKY effect arising in fourth order perturbation theory 
for two closely spaced dots and is thought to compete 
with Kondo physics. In particular, the RKKY-Kondo ef- 
fect is non-perturbative in virtual hopping processes. We 
find that the RKKY-Kondo effect is generically present 
provided edi 7^ ^d2 and so there are unequal numbers 
of electrons on each dot (as indicated in the cartoon in 
Figure 2). 

This antiferro-RKKY Kondo effect is, in a sense, un- 
surprising. The Friedel sum rule (FSR) dictates that the 
at the p.h.s. point, the scattering phase equal tt and 
G vanish. If we were instead to have a ferro-RKKY ef- 
fect, and so an underscreened spin 1 impurity, G would 
be maximal and so violate the FSR. Our finding of a 
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Fermi liquid fixed point is supported indirectly by earlier 
work on the two impurity Anderson model pij . Although 
there the focus is upon two electron channels coupled to 
two impurities, the finding is both that the physics is 
Fermi liquid and that the one to two channel crossover 
is smooth. We have verified singlet formation both by 
demonstrating from the Bethe ansatz equations that the 
entropy vanishes in the zero temperature limit and, sep- 
arately, from a slave boson mean field analysis. 

The Abrikosov-Suhl resonance associated with this 
RKKY-Kondo effect can be computed along the lines of 
Rcf. (up to a multiplicative constant): 

p(e) = cos(/37r)(r^^(e~2 + l))/(2e~2 cos(2/3^) + + 1), 

where e = e/TjiK and Trk is the RKKY Kondo tempera- 
ture, Tra' ^ VUi^i exp(— TrtZ/SFi). The parameter (3 is 
if Ae > ri,2 and is 1/2 if < Ae < Ti^2- As Ae 0, 
the resonance evolves from a Lorentzian centered at zero 
energy to a structure with split peaks. The scale Trk 
governs the leading corrections to the conductance as a 
function of Zceman field (H) and temperature (T). For 
Ae » ri^2, we compute these at the p.h.s. point to have 
the Fermi-liquid form (in units of 2e'^/h): G{T/Trk) = 
tt\T/Trk)V^ and G{H/Trk) = tt\H/Trk?/A. 

As we move away from the p.h.s. point through in- 
creasing the gate voltage, we begin to empty the dots. 
With the assumed asymmetry in e^i and 6^2, we arrive 
at a point where the dot system has roughly one electron 
sitting predominantly on dot 1. At this value of the gate 
voltage, we expect ordinary Kondo physics to be opera- 
tive. Defining Tk via its relation to the static impurity 
susceptibility, i.e. Ximp = (47V)~^, the corresponding 
leading contributions to the conductances (in units of 
2eV/i) are G(T/Tk) = 1 - {■K^/l%)^TKf and 
G{H/Tk) = 1 - (7rVl6)(ir/rif)2 [l^. From a numeri- 
cal analysis, we know that Tk has a single dot form, i.e. 
Tk ^ y!7T772exp(7redi(edi + U)/2UiTi) ^M,. 

As tdi is further increased, we see both a vanishing 
of the conductance and an unusual form of the FSR, a 
mark of the effects of interference. The vanishing of the 
conductance may reflect interference alone: it is present 
as well in the non-interacting case. However the form 
the FSR takes reflects both interference and interactions. 
The FSR relates the scattering phase to the number of 
displaced electrons, i.e. 5, 
be 



Ti'^disCT- '^discr IS defined to 



nda 



dx 



(10) 



and contains contributions from both the occupancy of 
the dots, Hda, and deviations in the lead electron density 
from coupling the dots to the leads 27[. As is evident 
in Figure 2, we have an unusual situation where both 
contributions to Jidiso-, and not merely Uda, are finite: 
Uda- is always manifestly positive while ridiscr is negative 
over a range of e^i . 



One last feature to the linear response conductance we 
wish to point out is the disappearance of the Abrikosov- 
Suhl resonance at precisely e^i = ed2- The transition is 
first order as Trk itself does not vanish as e^i — > ed2- 
The origin of this critical point lies in the decoupling 
of one dot degree of freedom (d.o.f) if e^i = ed2 (and 



3nly if), as can be seen via a change of basis d, 



e/o 



(Fi/adi ± T2/id2)/{rl + Vly'^. This discontinuous be- 
havior at edi = £d2 however can be transformed into a 
smooth crossover by weakly coupling a second channel of 
electrons to the dot (2^ and so recoupling the odd dot 
d.o.f. 

Apart from the behavior at edi = ^d2i we generally 
expect the above physics to be robust against small vi- 
olations of the integrability constraints (Eqn. 9) for a 
number of reasons [21|: i) the ground state of the dot- 
lead system is already robustly established (unlike when 
perturbation theory in V is done for a single dot-lead); 
and ii) a Schrieffer- Wolfe transformation in the Kondo 
regime is unaffected by (weak) violations of Eqn. (7). 
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